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Abstract. We present radiation-driven wind models for Luminous Blue Variables (LBVs) and predict their mass- 
loss rates. We study the effects of lower masses and modified abundances in comparison to the normal OB 
supergiants, and we find that the main difference in mass loss is due to the lower masses of LBVs. In addition, we 
find that the increase in helium abundance changes the mass-loss properties by small amounts (up to about 0.2 dex 
in log M), while CNO processing is relatively unimportant for the mass-loss rate. A comparison between our mass 
loss predictions and the observations is performed for four relatively well-studied LBVs. The comparison shows 
that (i) the winds of LBVs are driven by radiation pressure on spectral lines, (ii) the variable mass loss behaviour 
of LBVs during their S Doradus-type variation cycles is explained by changes in the line driving efficiency, notably 
due to the recombination/ionisation of Feiv/lll and Felll/ll, and finally, (iii) the winds of LBVs can be used to 
derive their masses, as exemplified by the case of AG Car, for which we derive a present-day mass of 35 Mq. 

Key words. Stars: early-type - Stars: mass-loss - Stars: supergiants - Stars: winds, outfiows - Stars: evolution 



1. Introduction 

Luminous Blue Variables (LBVs) are massive stars in a 
brief (~ 10^ yr) but violent post-main-sequence phase of 
evolution, during which they lose the remaining parts 
of their hydrogen-rich outer layers, before entering the 
Wolf-Rayet stage. They show a complex pattern of vi- 
sual brightness variations (Humphreys & Davidson 1994; 
Lamers 1995; van Genderen 2001), but most characteris- 
tic are their visual light variations of ~ 1 - 2 magni- 
tudes on timescales of years. These changes are referred 
to as "typical" or "S Doradus-type" outbursts. During 
these variability cycles, the stars move horizontally in the 
Hertzsprung-Russell Diagram (HRD) as they expand in 
radius at approximately constant luminosity. At minimum 
visual brightness the star is relatively hot (Teff ~ 20 000 - 
30000 K) and small, whereas the star is large and rela- 
tively cool (Toff ~ 8 000 K) at visual maximum. Although 
it has been established that the typical variations are re- 
lated to changes in stellar radius rather than the formation 
of a pseudo-photosphere (Leitherer et al. 1989; de Koter 
et al. 1996), the physical reason for these radius changes - 
which can be as large as factors of ten - is as yet uniden- 
tified (see Nota & Lamers 1997 for overviews). 

The strong stellar winds of LBV stars {M ~ 10"** - 
10~^ M0yr~^) show a wide variety of mass loss behaviour. 
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which is poorly understood. There are cases like R71 
where the mass loss increases while the star expands 
(Leitherer et al. 1989), whereas RllO shows a behaviour 
that is the exact opposite: while the star increases in size, 
its mass- loss rate drops (Stahl et al. 1990). These results 
imply that there is no straightforward correlation between 
mass loss and radius (or effective temperature) for LBVs 
(Leitherer 1997). 

Recent radiation-driven wind models for normal O and 
B stars however show that one does not expect such a 
strict monotonous behaviour of mass loss with effective 
temperature. Vink et al. (1999, 2000) have found an over- 
all trend that cooler stars - at constant luminosity - ex- 
perience smaller mass loss, but at specific effective tem- 
peratures the mass-loss rate jumps upward. The overall 
trend can be explained by a growing mismatch between 
the position of the bulk of the driving lines (mostly in 
the ultraviolet) and the location of the flux maximum, 
which shifts gradually towards the optical for cooler stars. 
Superimposed on this there are jumps, where M increases 
steeply due to recombinations of dominant line-driving 
ions, specifically those of iron. At Toff ~ 25 000 K, the mass 
loss is predicted to increase by about a factor of five due 
to the recombination of Feiv to Felii. Observational ev- 
idence for this so-called 'bi-stability' jump has been pre- 
sented by Lamers et al. (1995), who discovered an abrupt 
discontinuity in the carbon ionisation and the terminal 
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wind velocities of stars at spectral type Bl (Tcff ~21000 
K). They also present evidence for a second jump at spec- 
tral type AO (Toff ~ 10000 K), which is anticipated to be 
due to the recombination of Fe ill to Fe ii (Leitherer et al. 
1989). 

The main goal of this paper is to investigate whether 
radiation pressure on spectral lines is responsible for the 
mass-loss rates and mass-loss variability of LBV winds, or 
that alternative mechanisms, such as pulsations or mag- 
netic fields, need to be invoked. If radiation-pressure is 
the main driver, just as it is in OB-type stars, then one 
expects the effective temperatures at the epochs at which 
M has been deduced, to determine whether an increase or 
decrease in mass loss occurs. In such a picture, the seem- 
ingly contradictory mass loss behaviour of e.g. R71 and 
RllO may be explained in terms of their differing effec- 
tive temperature intervals with respect to the bi-stability 
jumps (de Koter 1997; Lamers 1997). 

The mass-loss rates of LBVs are expected to differ 
from those of normal OB stars - at the same position in 
the HRD - as the result of two effects. First, their stellar 
masses (and consequently their mass over luminosity ra- 
tios) are expected to be lower as LBVs are more evolved. 
Second, their surface abundances are expected to be differ- 
ent as products of nucleosynthesis may have reached the 
stellar surface. Observations show that LBVs are indeed 
enriched in helium and nitrogen, and depleted in carbon 
and oxygen (e.g. Smith et al. 1994; Najarro et al. 1997). 
The LBV mass-loss predictions presented in this paper 
take these effects into account. 

As LBVs have the property that their effective tem- 
peratures change substantially on timescales of years, they 
offer the unique opportunity to study mass loss as a func- 
tion of effective temperature for individual stars. In this 
respect it is interesting to note that Stahl et al. (2001) 
have recently observed and analysed the Galactic LBV 
AG Carinae and found an increase in the mass-loss rate 
by a factor of five during its excursion from minimum to 
maximum visual light. Intriguingly, this increment in mass 
loss is similar in magnitude to the predictions of OB star 
mass loss around the bi-stability jump (Vink et al. 1999). 

The paper is organised as follows. We first describe 
the Monte Carlo method used to predict mass-loss rates 
(Sect.H). In Sect.|^ we present the LBV mass-loss pre- 
dictions. A comparison between these theoretical mass- 
loss rates and all available mass-loss observations is sub- 
sequently presented in Sect. |. Special emphasis is at- 
tributed to AG Car (Sect.^^), which is by far the best 
studied member of its group. Interestingly, the set of Af- 
determinations for this star provides such tight constraints 
that, in addition, the stellar mass of AG Car can be deter- 
mined (Sect.^). We present our conclusions in Sect.||. 



2. Method to predict mass loss 



2.1. A brief description of the Monte Carlo method 

We use Monte Carlo calculations of photons travelling 
through an atmosphere that accounts for an outflowing 
wind to predict mass-loss rates. In essence, the mass loss 
follows from equating the gain of kinetic energy of the 
outflowing gas to the decrease of radiative energy. From 
an outside observer's point of view, this decrease of ra- 
diative energy occurs when photons transfer momentum 
(and energy) in interactions with moving ions in the flow. 
An iterative process is used to find the model in which 
the input mass loss is equal to that found from the global 
radiative to kinetic energy conversion. This M value is 
the predicted mass-loss rate. Details of the Monte-Carlo 
method are given in Abbott & Lucy (1985), de Koter et 
al. (1997) and Vink et al. (1999). The model atmospheres 
used in the process are calculated using the non-LTE code 
ISA-WIND and involve a gradual transition of density, con- 
sequently of velocity (through the equation of mass conti- 
nuity), and temperature between the stellar photosphere 
and wind region. The chemical species for which the statis- 
tical equilibrium equations are solved explicitly are H, He, 
C, N, O, and Si. For a full description of the code, we refer 
the reader to de Koter et al. (1993, 1997). The iron-group 
elements, which are important for the radiative driving, 
are not explicitly accounted for but are treated in a gen- 
eralised version of the "modified nebular approximation" 
described by Lucy (1987, 1999). 

Our method to predict mass loss is distinct from the 
approach in which the line force is parameterised in terms 
of force multiplier parameters (e.g. Castor, Abbott & 
Klein 1975, Pauldrach et al. 1994). The main advantage 
of our approach is that it accounts for multiple photon 
scattering processes, which have been found to be im- 
portant even for normal O- and B-type stars (e.g. Vink 
et al. 2000), but which are not treated in the line-force- 
parameterisation method. A disadvantage of our approach 
may be that we do not de facto solve the momentum equa- 
tion, i.e. the predicted mass loss values depend to some ex- 
tend on the adopted velocity stratification v{r). For this 
we assume a /9-type velocity law, which has been found to 
yield excellent results in modelling the spectra of hot star 
winds (e.g. Puis et al. 1996). We estimate the intrinsic ac- 
curacy of our Monte Carlo calculations to be ^ 0.05 dex 
in the predicted M. 



2.2. Assumptions in the models 

Modified nebular approximation. The simplified treatment 
of the ionisation and excitation state of the iron-group el- 
ements, i.e. by means of the modified nebular approxima- 
tion, may result in a (systematic) shift of the temperature 
at which the dominant ionisation state of these species 
changes relative to full non-LTE calculations. An indica- 
tion that such an offset may indeed be present has been 
presented by Vink et al. (1999). These authors identify the 
occurrence of strong changes in the terminal velocity and 
ionisation of the winds of OB supergiants at spectral type 
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Bl (i.e. at 21 000 K; Lamers ct al. 1995) to be a result of 
the transition of the dominant iron ionisation from Fe iv 
to III. Using a similar approach as the one used in this 
study, they predict a jump in the mass-loss rate as a re- 
sult of this recombination, but they predict the location of 
this jump to be at ~ 25 000 K, i.e. at an effective tempera- 
ture of ~ 4 000 K higher than the observed temperature of 
the jump. OB supergiant simulations over a wider range 
of wind densities indicate that the ionisation balance of 
Fe iv/iii is mostly sensitive to temperature, and not den- 
sity. Therefore, the bi-stability jump was found to occur 
over only a relatively narrow range of 22.5 < Tcff < 26 
kK (Vink et al. 2000). In other words, the offset between 
the predicted bi-stability jumps and the observed Bl jump 
was found to be in the range 1500-5 500 K, with 4000 K 
being a typical value. One may therefore expect that for 
the dense winds of LBVs, the offset will be at least of the 
same magnitude as the maximum B supergiant offset, i.e. 
5 500 K. We will later see (in Sect.^^ that the offset is 
~ 6 000 K. We note that for stars that vary their effec- 
tive temperatures in a range without iron recombinations, 
one does not expect the modified nebular approximation 
to affect the predicted mass loss. Indeed, in the case of 
0-type stars (where Fe iv is dominant over the full 0-star 
temperature range) good agreement was found between 
our predictions that use the modified nebular approxima- 
tions and the observed mass-loss rates (Vink et al. 2000; 
Benagha et al. 2001). 

In view of the above, the use of the modified nebular 
approximation is anticipated to cause an offset in the pre- 
dicted position of the bi-stability jumps. We will apply a 
corrective shift to the locations of the jumps in our mod- 
els to facilitate a meaningful comparison between these 
predictions and the observations (Sect.^). 

Stationarity. We assume a stationary stellar wind, i.e. 
both the mass-loss rate and the velocity stratification are 
time- independent. The assumption of stationarity may 
not be too bad an assumption for predicting the time- 
averaged mass-loss rates of O stars (see e.g. Owocki et al. 
1988). However, for LBVs this may not be fully correct 
as the timescale of photometric variability Tphot becomes 
comparable to, or shorter than, the typical How timescale 
of the wind. If we define the latter as the timescale Tdyn for 
a gas element to travel from the photosphere to 100 stellar 
radii - typically the extent of the Ha formation region in 
LBVs - we find at minimum visual brightness Tdyn ~ days, 
while at maximum visual light Tdyn ~ months. For LBVs, 
significant photometric changes may occur on timescales 
of Tphot ~ months, and therefore the assumption of a sta- 
tionary wind is first expected to break down for LBVs in 
their transition from maximum to minimum visual light. 
One should realise that such a breakdown would not only 
occur in wind models to predict mass loss, but it would 
also affect the analysis of LBV spectral observations, e.g. 
the modelling of the Ha line profile in order to determine 
the mass-loss rate. A stationary model atmosphere cannot 



Table 1. Adopted stellar parameters for the grid of LBV 
models presented in this study. 



logL* 


M, 


To 






(Mq) 




(kK) 


5.5 


30 - 10 


0.18 - 0.53 


10.5 - 30.0 


6.0 


35 - 25 


0.48 - 0.67 


10.2 - 30.0 



predict the line cores (which are formed at low projected 
velocities) and the line wings (formed at high velocities) 
at the same time, as this would introduce significant errors 
in the deduced M values. 

For these reasons we focus our comparison of predic- 
tions with observations on the transition of LBVs from 
minimum to maximum visual li ght, and we do not include 
the reverse trajectory (see Sect. 1.2). 



Sphericity. Spectropolarimetric observations of the LBV 
AG Car and other LBVs have shown that the winds of 
these stars are not spherically symmetric (see e.g. Schulte- 
Ladbeck et al. 1994). This may have implications for both 
our mass-loss predictions, as well as the modelling of the 
Ha lines to deduce the observed mass-loss rates. As there 
is little quantitative information on the density contrast 
between the pole and equator in LBV winds, it remains 
to be seen how important these non-sphericity effects are 
with respect to the global wind properties. In any case 
there is no particular reason why a possible systematic 
non-sphericity effect would appear more prominent in one 
visual brightness state than in the other. 

2.3. The adopted parameters 

Luminosities and IVlasses. We calculate mass-loss rates 
for two series of LBV models, one having a luminosity 
log(L/LQ) = 5.5 and the other 6.0. The luminosities of 
most LBVs fall in this range, e.g. R71 and RllO are 
representative for the low luminosity value, while AG Car 
and R127 are representative for the higher one. Since the 
masses of LBVs are poorly known, we adopt a range of 
masses for both series (see Table 0). For the series with 
log(L/LQ) — 5.5, mass loss was computed for M — 30, 
25, 20, 15, 12 and 10 Mq. For the series with log(i/LQ) 
= 6.0, the adopted masses were M — 35, 30, and 25 Mq. 

The luminosity-to-mass ratio (L/M) can be expressed 
in terms of the Eddington factor Fc. This factor is defined 
as the ratio between the Newtonian acceleration and the 
radiative acceleration due to electron scattering, and is 
given by 



incGM 



7.66 10^'Ve 



L 



M 



(1) 



where the constants have their usual meaning and Ue is the 
electron scattering cross-section per unit mass. The value 
for (To depends on temperature and chemical composition. 
The range in Fc relevant for our model grid is given in 
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Table 2. Solar versus adopted LBV abundances. The lat- 
ter show surface helium enrichment and CNO processed 
material. 



Element 


Solar 


LBV 




(mass fraction) 


(mass fraction) 


H 


0.68 


0.38 


He 


0.30 


0.60 


C 


2.9 10"^ 


2.9 10"" 


N 


9.510"'' 


9.510"^ 





7.710"^ 


7.710"" 



Col. 3 of Table |l|. Note that Fc is considerably smaller 
than unity, signifying that the classical Eddington limit is 
not exceeded. The reason for this relatively low Fe is that 
although the masses are lower than for normal OB super- 
giants, the value for the electron scattering cross-section 
(To is also lower. This is a result of helium enrichment in 
LBV atmospheres. 

Effective Temperatures. Models have input temperatures 
computed for effective temperatures between 11000 and 
30000 K with a stepsize of 2 500 K in the range 12 500 
- 30 000 K. This input temperature Tin is defined by the 
relation: 

L = 4^i?f„allt (2) 

where i?in is the inner boundary of the atmospheric model. 
i?in is typically located at Rosseland optical depth 20 ^ 
tr ^ 25. Note that the actual effective temperature Toff, 
which follows from the model computation, will be some- 
what lower. The effective temperature is defined at the 
point where the thermalisation depth measured in the cen- 
ter of the V band (at 5555 A) equals l/\/3 (see Schmutz 
et al. (1990) and de Koter et al. (1996) for a detailed dis- 
cussion). For stars with relatively low mass loss and con- 
sequently mostly optically thin winds, such as the winds 
from normal OB supergiants, there is no difference be- 
tween the input Tin and Toff. But in the case of LBVs, 
which may lose mass at rates of up to M ~ lO"''M0yr"^, 
there is an offset between the two temperature values, but 
this offset is relatively small. Note that the discrepancy be- 
tween Tin and Toff is more dramatic for Wolf-Rayet stars, 
which have much denser winds. 

Abundances. Both spectroscopic analyses and models for 
stellar evolution show LBVs to be enriched in helium 
(Smith et al. 1994; Meynet et al. 1994; Najarro et al. 
1997). In addition to helium enrichment, CNO processed 
material that has reached the surface causes an increase in 
the amount of nitrogen. This occurs at the cost of oxygen 
and carbon, which become depleted. The total amount of 
CNO remains nevertheless constant. The relative changes 
in the elements CNO have been adopted from the abun- 
dance computations for massive stars with rotation in- 
duced mixing by Lamers et al. (2001) and have been ap- 



plied to the solar CNO abundances given by Anders & 
Grevesse (1989). Table ^ presents both the initial solar 
abundances (Col. 2) as well as the abundances modified 
by the CNO process (Col. 3). The metal abundance Z 
is kept constant at the solar value, for models computed 
with all species included. 

Velocity Laws. We have calculated M for wind models 
with a /3-type velocity law for the supersonic part of the 
wind: 

v{r)^voo{l-^^^ (3) 

where Voo is the terminal flow velocity. Below the sonic 
point, a smooth transition from this velocity structure is 
made to one that corresponds to the photospheric density 
structure and adopted M (see de Koter et al. 1997 and 
Vink et al. 1999 for details). A value oi [3 — \ was assumed 
in all calculations. The predicted mass-loss rates are found 
to be insensitive to the wind acceleration parameter when 
(3 is in the range 0.7 - 1.5 (Vink et al. 2000). We specified 
the terminal velocity adopting for the ratio Voa/vesc the 
values 1.3, 2.0, and 2.6. Lamers et al. (1995) found that a 
ratio of 2.6 for Galactic supergiants of spectral type earlier 
than Bl; a ratio of 1.3 for Bl to AO stars, and a tentative 
ratio of 0.7 for later spectral types. 

3. The predictions of mass loss of LBVs 

Differences in the wind properties of an LBV star and a 
normal O- or B-type supergiant in the same location of 
the HRD are expected to be the result of differences in 
stellar mass (specifically in L/M) and in surface chemical 
composition. Regarding the latter, in LBVs helium may be 
enriched at the surface and CNO-cycle processed material 
may have appeared. This yields an increase in N and a 
decrease in C and O relative to the initial abundances. We 
will now discuss the difference in mass-loss rate between 
LBVs and OB-stars as a result of these effects. 

3.1. The effect of CNO processed material on M 

Figure |l| shows a comparison in mass-loss behaviour be- 
tween models with initial (solid lines) and CNO processed 
abundances (dashed lines). The results are very similar, 
i.e. the occurrence of CNO processed material does not 
change the predicted mass loss in any major way. There 
are two reasons for this. First, in the case of Galactic 
abundances (which are being discussed here) the domi- 
nant wind driving element is iron, relative to which CNO 
only plays a modest role in determining the mass-loss rate 
(Vink et al. 1999, Puis et al. 2000). Second, as the CNO 
elements only serve as a catalyst in the nuclear reactions, 
the total number C-|-N-f O remains constant. The net ef- 
fect of the nitrogen enrichment and the carbon/oxygen 
depletion on the line driving is therefore expected to be 
small, as it relies on differential effects in the number of 
lines provided by the relevant ions of these species. 
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Fig. 1. Predictions of mass loss as a function of effective temperature for a star ■with solar (dashed lines) versus CNO 
processed surface abundances (solid lines). As CNO serves as a catalyst in the nuclear reactions the total number 
C+N+0 atoms is invariant. The adopted stellar parameters are log(L/Lo)= 6.0; M = SOM©; (X,Y) = (0.68,0.30). 
Adopted values for the ratio bet'ween terminal and effective escape velocity are given in the legend. 



3.2. The effect of helium enrichment on M 

To investigate the difference in mass-loss behaviour as 
a result of changes in the hydrogen and helium abun- 
dance, we compared predictions for (X,Y) = (0.68, 0.30) 
and (0.38,0.60). The models have log(L/Lo)= 6.0 and M 
= 60 Mq, and are presented in Figure |[ The solar abun- 
dance models (dashed lines) have been taken from Vink et 
al. (2000). The models -with enhanced helium (solid lines) 
are characteristic for ■what is expected in the LBV phase 
(Meynet et al. 1994). 

The most pronounced difference between the two sets 
of models is that at T ^ 25 000 K (for Woo/wosc = 1-3 and 
2.0) the hehum enriched models show a higher mass loss, 
while at lower temperatures they feature a lower mass 
loss. The difference in M extends up to ~0.2 dex. A sec- 
ond difference is that the position of the bi-stability jump 
at ~ 25 000 K (where the iron ionisation in the region of 
strong wind acceleration flips between Feiv and ill) ap- 
pears to shift slightly towards higher temperature, while 
the opposite occurs for the jump at ^ 15 000 K (i.e. where 
Felii changes over to Fell). 

The cause of this behaviour is not the result of a differ- 
ence in the line driving of H i relative to He I and He ii, as 
hydrogen and helium are both irrelevant for the line driv- 



ing force at Galactic abundances (Vink et al. 2001). As 
it is the iron lines that dominate the wind driving (below 
the sonic point), the relevant question concerns the domi- 
nant ionisation stage of iron, and how well the wavelength 
positions of the Fe lines match the spectral region where 
the bulk of the flux is emitted. The many lines of Feiv 
and Fe ill are however spread throughout the Lyman and 
Balmer continua, and do not show a particularly strong 
preference for either side of the Lyman jump (at 912 A). 
Therefore, the cause of the differences in mass loss be- 
haviour must instead be related to changes in the under- 
lying flux distribution, in particular the strength of the 
Lyman jump. 

When hydrogen is (partly) replaced by helium, the 
size of the Lyman jump decreases. This implies that more 
photons are available in the Lyman continuum (at wave- 
lengths below 912 A). Since the luminosity of the models 
is kept constant, less flux will be emitted in the Balmer 
continuum. At relatively low Teff, the bulk of the flux is 
emitted in the Balmer continuum, therefore one expects 
the mass loss in the helium-enriched models to be lower; 
at relatively high Toff, most of the flux is emitted in the 
Lyman continuum, and as the flux in the helium-cnrichcd 
models is higher here, more mass loss occurs. 
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Fig. 2. Predictions of mass loss as a function of effective temperature for a star with solar (dashed lines) versus 
helium enriched (X=0.38, Y=0.60) abundances (solid lines). The adopted stellar parameters are log(L/L0) = 6.0 and 
M = 60Mq. Adopted values for the ratio bet-ween terminal and effective escape velocity are given in the legend. 



The small differences in the Toff positions of the t-wo 
bi-stability jumps can be explained in a similar -way. In 
the regime of the Feiv to ill jump at ^25 000 K, the 
wind densities in the helium-enriched models are higher, 
therefore recombination/ionisation will occur at a slightly 
higher temperature relative to the solar abundance mod- 
els. For the Feiii/ii jump at ~ 15 000 K the situation is 
reversed: the wind densities in the helium-enriched mod- 
els are lower, therefore the jump occurs at somewhat lower 
Tcs- 

3.3. The effect of lower stellar masses on M 

The effect of decreasing the stellar mass is shown in Fig. || 
for a set of models that have log(L/L0) = 5.5. The cal- 
culations show that for all temperatures the mass loss in- 
creases with decreasing stellar mass. For this luminosity, 
Vink et al. (2000) derived logAf cx — 1.3 logAf for masses 
in the range between 30 and 50 Mq. Here, we find that 
the predicted dependence is given by 



logM cx -1.81 (±0.05) logM 



for masses from 30 down to 10 Mq. For this mass range, 
this steep dependence implies an increase of the mass loss 
by a factor of seven, when the mass is decreased by a factor 



of three (from 30 to 10 Mq). This strong sensitivity can, 
in principle, be used as a diagnostic to derive current LBV 
masses. We will return to this in Sect. 4.2. 

Also visible in Fig. || is a gradual shift of the location 
of the bi-stability jumps towards higher temperature for 
lower stellar masses. This is as one would expect: the mod- 
els with low stellar mass have a relatively high mass loss, 
and consequently high wind density. The Saha-Boltzmann 
ionisation equilibrium (as formulated in the modified neb- 
ular approximation) then implies that recombination will 
occur at a somewhat higher temperature. 

For the series with luminosity log(L/LQ) — 6.0 pre- 
dictions for the mass-rate as a function of stellar mass are 
presented in Fig. ^. The calculations for 60 and 80 Mq 
shown in Vink et al. (2000) have been extended here to 35 
(top panel), 30 (middle panel), and 25 Mq (lower panel). 
Results are shown for three different ratios of Wqo / Vcsc- The 
same trends are visible and the same mass-loss vs. stellar 
mass dependence can be inferred. 



(4) 3.4. Summary: 



We conclude from the comparison between the mass loss 
of normal OB supergiants and LBVs that: (1) CNO pro- 
cessed material at the surface of LBVs does not result in 
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Fig. 3. Predictions of mass loss as a function of stellar mass for a star ■with solar abundances. The adopted luminosity 
is log(L/L0) — 5.5. We find that log A/ cx —1.8 logil/. The different masses in Mq are given in the plot. 



a noticeable change in M as the total number of effective 
driving lines is not significantly affected. (2) Helium en- 
richment in LBVs causes a modest change in mass loss 
■which can be traced back to a decreased Lyman jump at 
912 A. (3) The reduced stellar mass of LBVs causes a 
strong increase in the mass-loss rate relative to normal 
OB supergiants. 

4. Comparison between LBV predictions and 
observations 

4.1. Observed mass-loss behaviour of LBVs 

In a recent census, van Genderen (2001) identifies a to- 
tal number of 34 LBV members in the Galaxy and the 
Magellanic Clouds. He divides the LBVs into three groups: 
strong-active members, characterised by light amplitudes 
of AV > 0.5 mag; weak-active members sho'wing AV < 
0.5 mag; and ex- and dormant members ■which have not 
sho'wn any typical photometric variability in the previous 
century. As our interest is in the mass loss behaviour of 
LBVs as a function of effective temperature, ■we are espe- 
cially interested in the strong-active members. Although 
all of these 14 stars have been monitored photometri- 
cally over the last decades, only a handful have been sub- 
ject to quantitative spectroscopic investigation at different 
epochs of photometric variability. This implies that esti- 



mates of the basic stellar parameters, such as the mass- 
loss rate, are rarely found at different epochs. In this re- 
spect, AG Car is by far the best investigated member, 
•with detailed studies undertaken by Leitherer et al. (1994), 
Schmutz (1997) and most recently Stahl et al. (2001). We 
will focus on AG Car in Sect. 4.2, but we first present a 



brief overvie^w of the other three LBVs that have been 
investigated spectroscopically. These LBVs are all in the 
Large Magellanic Cloud. 



R71: Spectra of Radcliffe 71 at minimum and maximum 
visual light have been discussed in detail by Wolf et al. 
(1981). They derive temperatures for the t-wo phases of 
Teff ~ 13 600 and 6 000 K respectively. The luminosity of 
the star is about 10^'^ L©. From a simple H/3 line profile 
analysis they determine M — 10^^'^ Moyr~^ for the hot 
phase. For the cool phase they estimate from the observed 
infrared excess that Al = lO"'*''^ M0yr~^. In a historical 
context, it is essentially this difference of a factor 100 - 200 
that spurred the idea of explaining the radius changes of 
LBVs by a variable mass loss, leading to the formation of 
pseudo-photospheres at maximum light (Davidson 1987, 
Lamers 1987). However, for two reasons we suspect that 
the derived maximum mass loss may be incorrect. First, 
dust emission contributes to the near-infrared flux of R71 
(e.g. Voors et al. 1999) which Wolf et al. (1981) have not 
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loss of a factor 100-200. In case H/3 is optically thin, 
holds that 



it 



1/2 



(5) 



where T'^ denotes the continuum flux at 4861 A and W^'^ 
represents the equivalent width of H/3. The above equation 
assumes that the terminal velocity does not vary between 
minimum and maximum phase, which has indeed been 
found to be the case (Wolf et al. 1981). To estimate the 
mass loss ratio between visual minimum (min) and maxi- 
mum (max) light, we use the H/3 peak strength ratios, as 
given by Leitherer et al. (1989), of respectively 1.3 and 5.0 
to approximate the equivalent width ratio. We estimate 
the continuum flux ratio using the total visual magnitude 
change of Ay ~ 1 mag. We then arrive at a mass loss at 
maximum visual light of ~ 3 Mmin . In case the line would 
be optically thick, and the power-law index in Eq. (^ is 
alternatively given by 4/3 (de Koter et al. 1996), we find 
that Mniax ~5.5Mmin. Notc that model calculations of 
Leitherer et al. (1989) also indicate that the change in the 
H/3 profile points to a much more modest difference in 
mass loss between the two phases than the factor 100-200 
estimated by Wolf et al. (1981). We come to the conclusion 
that R 71 shows an increase in mass loss of a factor of 3 to 
5 when it changes its temperature from ^ 13 600 to 6 000 
K. Our models indeed predict a change of this magnitude 
as the star crosses the Feiii to Fell recombination jump 
at 10 000 K. 
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Fig. 4. Predicted mass-loss rates as a function of effective 
temperature for three adopted stellar masses. All models 
have log (L/L0) = 6.0. The values for the ratio between 
terminal and effective escape velocity are given in the up- 
per panel. 



corrected for. This implies that they have overestimated 
M . Second, the change in the H/3 profile between minimum 
and maximum light does not corroborate a change in mass 



R 127; R127 has been extensively monitored by the 
Heidelberg group (Stahl et al. 1983, Stahl & Wolf 1986, 
Wolf et al. 1988). Based on spectral type determination, 
the stellar temperature is found to vary between about 
30 000 and 10 000 K. The authors note a rather capri- 
cious behaviour of the mass-loss rate. We anticipate such 
behaviour as two iron recombination jumps are located 
within this temperature range, i.e. Feiv to ill at ~ 21 000 
K and Feiii to 11 at ~ 10 000 K. Unfortunately, a lack of 
well determined mass-loss rates at different phases in its 
variability cycle prevents a quantitative comparison. 

R 110: The lightcurve of R 110 has been discussed in detail 
by van Genderen et al. (1997). The star has been inves- 
tigated spectroscopically at two different epochs by Stahl 
et al. (1990). They found that in the course of brighten- 
ing, when the star changed its spectrum from B-type with 
P Cygni lines to normal FO, the mass loss decreased from 
an estimated few times 10~^ Moyr"^ to a value that can- 
not be much higher than about 10~^ Moyr"-'^. The impli- 
cation of the finding that M is smaller for lower Tes is that 
the cause of the excursions of LBVs in the HR-diagram 
cannot be the formation of an opaque stellar wind as the 
result of an increased mass loss. Our mass-loss predic- 
tions do indeed show Teff intervals for which the mass loss 
increases inversely proportional to the temperature. The 
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Fig. 5. Mass- loss rates of AG Car as a function of effective 
temperature as derived from analysis of the Ha profile by 
Stahl et al. (2001). The solid Hne connects points on a 
rising branch in the lightcurve, i.e. when the star is on its 
way from minimum to maximum visual light, and covers 
the epoch Dec 1990 - Feb 1995. The dashed line indicates 
the mass-loss behaviour on the falling branch, covering the 
period Feb 1995 ~ Aug 1999. 



uncertainty in the derived mass loss by Stahl et al. (1990) 
is however too large to allow for a meaningful quantitative 
comparison. 

4.2. The mass-loss behaviour of AG Car 

AG Car is by far the best observed LBV in terms of its 
mass-loss behaviour. Leitherer et al. (1994) derive M{Tcs) 
for the rising branch in its lightcurve covering the period 
from December 1990 to June 1992. Stahl et al. (2001) also 
investigate the mass loss behaviour of the star covering the 
rising and falling branches in the epoch December 1990 
- August 1999. The methods in these papers to deter- 
mine M are to some extent rather similar. In both stud- 
ies the luminosity of the star is assumed to be constant 
at log(L/L0) = 6.0, and the temperature is derived by 
fitting the observed flux in the 1^-band to that of non- 
LTE model atmospheres (see the respective papers for de- 
tails). The two methods differ in the way the mass loss 
is derived. Leitherer et al. used the equivalent width of 
Ha, while Stahl et al. modelled the Ha profile at selected 
phases. The latter method is more accurate, and it yields 
more reliable M values; it additionally allows one to set 
constraints on the velocity law. We will therefore com- 
pare our predictions to the data of Stahl et al. (2001). 
Their derived mass-loss rates are plotted against effective 
temperature in Fig. ^ for both the rising (solid line) and 
the falling branch (dotted line). Note that there is a dif- 
ference in behaviour between the star crossing over from 
minimum to maximum visual light, and visa versa. Our 
models do not predict such behaviour, as the basic stellar 
parameters (including stellar mass) are identical at similar 



positions in the HRD on the rising and the falling branch. 
We suspect that the differences in M between the two 
branches, as derived by Stahl et al. (2001), are the result 
of the breakdown of the assumption of a stationary wind 
(see Sect. 2.2) for the falling branch. The long dynamical 
timescale of the flow at maximum light implies that mate- 
rial lost in that phase may still be in close vicinity to the 
star, affecting the mass loss determination in the falling 
branch leading to a spuriously high (surface) mass-loss 
rate. A second explanation may be the release of gravi- 
tational energy when the star returns from maximum to 
minimum. If this plays a role, the assumption of constant 
bolometric magnitude during the variability is no longer 
valid (see Lamers et al. 1995). Whether or not AG Car re- 
mains at constant luminosity during its typical variations 
has however not yet been firmly established. Because of 
the two above-mentioned complications, we will focus the 
comparison of our predictions to the observed mass-loss 
rates to the rising branch only. 

In Fig. 1^ we compare our mass-loss predictions for 
AG Car with the rates determined by Stahl et al. (2001) 
for the rising branch over the epoch Dec 1990 - Aug 1999, 
where the temperature decreases from 24 000 to 9 000 K. 
In our calculations we adopt log(L/L0) — 6.0; a stellar 
mass M = 35 M©; a helium mass fraction Y = 0.60, and 
a constant ratio of terminal over escape velocity Woo/'^esc 
— 1.3. The luminosity and helium abundance are very 
similar to that assumed by Stahl et al. (2001). The stellar 
mass is a free parameter or, phrased differently, fitting the 
mass loss allows one to determine the mass (see Sect. 4.3). 
Unfortunately, the ratio Woo/^osc turns out to be rather 
poorly constrained by observations, as Ha only allows to 
determine a lower limit to Voo (Stahl et al. 2001). Taking 
these lower limits and correcting the escape velocity for 
electron scattering, Stahl et al. find lower limits in the 
range 0.2 ^ Woo/wesc ^1.1- We have opted to use our low- 
est calculated Woo/^^csc for the whole temperature range, 
as this is the most appropriate value. We need to make 
two remarks with respect to this choice. 

First, the adopted ratio of the terminal over the effec- 
tive escape velocity influences the present-day mass de- 
termination of the star. For instance, if we would adopt a 
value Vca/vesc = 2.0 we would find M = 30 Mq, while a 
choice of Voo/vesc = 2.6 would yield M = 25 Mq. To first 
order, this positive relation between the velocity ratio and 
the stellar mass is explained as follows. For fixed stellar 
parameters, the photon momentum (L/c) transfered to 
the wind (Mvoo) is more or less constant. It would there- 
fore be necessary to decrease the stellar mass to match 
the observed mass-loss rates. However, if the actual value 
of Voo/vcsc were less than 1.3, we would anticipate that 
the mass would be lower than 35 Mq rather than higher. 
This because Doppler shifting of lines is strongly reduced 
for low velocity winds, which decreases the effectiveness 
of wind driving and would indeed require a lower stellar 
mass to match the observed M. 

Second, the temperature range investigated by Stahl et 
al. (2001) includes the regime where iron recombines from 



10 



Jorick S. Vink & Alex de Koter: Predictions of variable mass loss for LBVs 



© 



o 



■3.5 



■4.0 



-4.5 



■5.0 




30 



25 



20 15 
Teff (kK) 



10 



5 



Fig. 6. Comparison of the predicted (dotted line) and observed (dashed line) mass loss of AG Car. The adopted 
stellar parameters are \og{L/LQ)= 6.0; (X,Y) = (0.28,0.60), and Uoo/wosc = 1-3. All calculations have been shifted by 
ATeff = —6 000 K to correct for a systematic offset in the predicted ionisation balance of iron. The good agreement 
allows us to determine the present-day mass of AG Car to be 35 Mq by fitting the mass-loss rates. 



Felll to Fell. As a result of this, we predict a jump in mass 
loss around 15 000 K. Possibly, this jump is related to the 
drop in the ratio Woo/^^csc from about 1.3 to 0.7 at spectral 
type AO, as identified by Lamers et al. (1995) on the ba- 
sis of observed terminal velocities of supergiants. Because 
AO I stars have TcS ~ 10 000 K, a similar shift bet'ween pre- 
dicted and observed ionisation appears to be present as 
■was found for the recombination of Fe I'V to Fe ill around 



spectral type Bl (cf. Sect. 2.2). In order to compensate for 
this offset, we shift our calculations by ATcff = — 6 000 K. 
The observed change in terminal velocity at spectral type 
AO -would affect the predicted M for the points at visual 
maximum light (i.e. at Toff ^10 000 K). As this is at the 
very lo^w end of the investigated temperature range, we 
have opted not to decrease Woo/fosc further, i.e. below 1.3. 

Figure ^ shows that after accounting for the corrective 
shift AToff, the observed and predicted mass loss agree 
within about 0.1 dex. Taking into account that M{Tcff) 
shows a complex behaviour, fluctuating over more than 
0.5 dex, this is a surprisingly good result. It shows that the 
mass loss variability of AG Car is the result of changes in 
the ionisation of the dominant line-driving element, iron. 



4.3. The stellar mass of AG Car 

The M{Tcs) behaviour can also be used to constrain the 
stellar mass. The reason is that the predicted mass-loss 
rates are not only a function of T^s, but also of stel- 
lar mass (see Eq. ^) . Comparing the AG Car data with 
M(Teff) predictions for different masses results in a best 
fit for the present-day mass of 35 ± 5 Mq . Note that this 
mass estimate does not account for possible systematic er- 
rors. Notably, one should be aware that the Ha spectral 
modelfing procedure of Stahl et al. (2001) did not account 
for clumping. In the case of Wolf-Rayet stars it is well- 
established that the neglect of a clumping correction may 
lead to systematic errors in the mass-loss rates by factors 
of ~ 2-4 (see e.g. Hillier 2000). If clumping also plays 
an important role in the winds of LBVs, this would im- 
ply that the mass-loss rates derived by Stahl et al. are 
too high. Consequently, our mass derivation would be un- 
derestimated. We will now discuss whether a present-day 
mass of 35 M0 is a reasonable value for AG Car. 

Applying multiple distance indicators, Humphreys et 
al. (1989) constrain the distance to AG Car to be 5 - 7 
kpc, implying a luminosity of log(L/L0) — 6.22 ± 0.16. 
Comparison with evolutionary tracks (Schaller et al. 1992) 
indicates that the most appropriate initial mass for such 
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a large luminosity equals 120 Mq (but it would be only 
85 Mq when main sequence mass loss is taken twice the 
canonical value; Meynet et al. 1994). The extreme bright- 
ness of such a star would imply a very high mass loss, 
which affects the evolution in a profound way. The evolu- 
tionary tracks show that a star with an initial mass of 120 
Mq evolves to the Wolf-Rayet phase almost directly from 
the main sequence, without even entering a post-main se- 
quence supergiant phase. This scenario is very difficult 
to reconcile with the (range of) spectroscopic properties 
observed for AG Car, which indicates that its luminosity 
must be less than log(L/L0) = 6.22. For these reasons, 
we adopt \og{L/LQ) = 6.0, a value which is still in fair 
agreement with the lower bound found by Humphreys et 
al. (1989). The corollary of all this is that the initial mass 
for AG Car is only 65 Mq, if we assume that the star is 
indeed in a post-main-sequence phase. 

Voors ct al. (2000) have determined the mass of the 
dust shell around AG Car and arrived at 0.25 Mq, with an 
accuracy of about a factor of two. This error is dominated 
by the uncertainty in the maximum size of the dust grains. 
Adopting the usual gas/dust ratio of 100, one finds a total 
circumstellar mass of about 25 Mq. If we assume that 
this material is of stellar origin, which is supported by 
the enhanced N abundance found in the ionised part of 
the nebula (Hutsemekers & van Drom 1991; Smith et al. 
1998), the upper limit to the present-day mass of the star 
is '~4OM0. This is in good agreement with the mass of 
35 Mq that we have derived in this study. 

5. Summary and Conclusions 

We have presented radiation-driven wind models for a set 
of stellar masses and luminosities characteristic for LBVs 
and have predicted their mass- loss rates. We studied the 
effects of lower masses and modified abundances on the 
mass-loss rates of LBVs in comparison to normal OB su- 
pergiants. We have found that the main difference is due 
to the lower masses of LBVs. In addition, we have found 
that an increased helium abundance changes the mass loss 
properties by small amounts (of up to about 0.2 dex in log 
M), while CNO processing is relatively unimportant for 
the mass-loss rate. 

The comparison of our predictions with the observed 
mass- loss behaviour of LBVs has shown that: 

1. The dominant driving mechanism of LBV winds is 

radiation pressure on spectral lines. Other mechanisms 
are expected to play only a minor role. 

2. The observed variable mass loss behaviour of LBVs of 
up to ~ 0.5 dex during their S Dor adus- type variation 
cycles can be explained in terms of changes in the 
efficiency of line driving, as a result of the recombina- 
tion/ionisation of Feiv/iil and Feiii/il. 

3. The M (Tcff) behaviour of AG Car can be matched 
to within 0.1 dex when a shift in temperature of 



ATeff = —6 000 K is applied to our predictions. This 
shift is required to correct for an inaccurate calcula- 
tion of the ionisation balance of iron, as we have used 
a modified nebular approximation for this element. 
Note that the shift is consistent with constraints set 
by observations of supergiants which indicate that the 
Feiv to III and the Feiii to ii recombinations indeed 
occur at spectral types Bl and AO (Lamers et al. 1995). 

4. We have deduced a value of 35 Mq for the present-day 
mass of AG Car. If we adopt a luminosity of 10^ 
Z/0 (which implies an initial mass of ~65M0) and 
correct for the amount of ejected material found in 
the circumstellar environment of the star (~25M0; 
Voors et al. 2000), we arrive at a present-day mass 
upper limit of ~4OM0, which is consistent with our 
finding of 35 M0. 
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